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We analyze the performance of quantum ratchets by considering the dynamics of an initially 
localized wave packet loaded into a flashing periodic potential. The directed center-of-mass motion 
can be initiated by the uniform modulation of the potential height, provided that the modulation 
protocol breaks all relevant time- and spatial reflection symmetries. A poor performance of quantum 
ratchet transport is characterized by a slow net motion and a fast diffusive spreading of the wave 
packet, while the desirable optimal performance is the contrary. By invoking a quantum analog of 
the classical Peclet number, namely the quotient of the group velocity and the dispersion of the 
propagating wave packet, we calibrate the transport properties of flashing quantum ratchets and 
discuss the mechanisms that yield low-dispersive directed transport. 
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I. INTRODUCTION 

Long-range, lasting and controllable quantum trans- 
port in optical and magnetic potentials is a key to all- 
atom optical and atom-on-chip devices. The ratchet ef- 
fect [il-Q presents a suitable toolbox for manipulations 
of ultracold matter, which allows to set a quantum par- 
ticle into directed motion without applying gradients 
or running- wave potentials [13] ■ Periodic modulations 
of the confining potential is prerequisite for the effect 
to occur, and there are plenty of different setups and 
blueprints of ratchet machinery [l], [ll|. An interest- 
ing class of ratchet systems are dissipation-free Hamil- 
tonian ratchets, classical Id Il2| and, even more in- 
triguingly, quantum ones [l|, lZHSI, Il3l4l8| . For exam- 
ple, a bi-harmonic flashing potential with periodically 
and uniformly modulated potential height |18j . could be 
turned into a conveyer belt for a Bosc-Einstcin conden- 
sate (BEC) of rubidium atoms fioj . 

The quantum ratchet effect can be utilized as a promis- 
ing device capable to deliver ultracold atoms to desirable 
locations. Yet the coherent quantum transport in ac- 
driven periodic potentials is severely hampered by diffu- 
sion, additionally enhanced by tunneling effects [23| • The 
issue of the transport efficiency [2l| - |23| becomes now of 
importance, meaning that one should search for a set of 
optimal parameters in order to maximize the correspond- 
ingly chosen efficiency measure. It is intuitive that the 
ratchet transport with large transport velocity and mini- 
mal dispersion rate would be preferable in the context of 
the delivery problem. This reasoning naturally leads to 
a q uantum analog of the classical Peclet number concept 
[2J], which has already been used in the field of classi- 
cal ratchets, both for overdamped [1^ and, as well, 
for underdamped [27| models. However, the extension 
of this classical measure to the quantum limit is by no 
means straightforward. Different from the classical limit, 
where the diffusive dispersion scales linearly, (t^(<) ~ t, 
the quantum diffusion is ballistic-like, cr^(i) ~ (20| . 
This fact makes the direct implication of the Peclet num- 
ber for the quantum ratchet problem not feasible. In 



this work we introduce a quantum analog of the Peclet 
number in order to specify the quality of ratchet trans- 
port. Using this concept we then aim in exploring the 
impact of diffusion and tunneling on the nonequilibrium 
quantum transport, focusing on initial conditions in the 
form of a localized Gaussian-like wave packet. Our ob- 
jective is to identify optimal regimes, when the directed 
transport is minimally swamped by the dispersive spread- 
ing. By use of the Floquct formalism [28j , which provides 
the complete information on the state of the system at 
any instant of time, we demonstrate how the analysis of 
the Floquet spectrum of the driven system allows one 
to judge the quality of a flashing ratchet as a suitable 
delivery vehicle for ultracold atoms. 



II. HAMILTONIAN QUANTUM RATCHET 
SETUP 

We consider a quantum particle loaded into a periodic 
potential, which is periodically modulated in time [l8| . 
The dynamics of the system is governed by a Schrodinger 
equation. 



ih^^m)) = H{t)m)), 



where the Hamiltonian Hit) reads 



H{x,p,t) 



P 

2m 



U{x)E(t). 



(1) 



(2) 



The flashing potential is formed by a bi-harmonic lattice 
of spatial period L, U{x) = U{x + L), 



U{x) = KU(3[cos{kx) s cos{2kx + 6p% k 



27r 



(3) 



whose amplitude is modulated by a bi-harmonic driving 
function, 



E{t)=E{t-to) = 
El cos[w(t - to)] + E2 cos[2a;(t - to) + S] 



(4) 
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symmetry 


transformation 


validity condition 


Si 


{x,p,t,K) {-X,-p,t,-K) 


U{x) = U{-x) 
requires 

s = or Op — ztniTT, m = 0, 1, ... 


52 


{x,p,t,K) -> {x,-p, -t,-K) 


E{t) = E{~t) 
requires 

E2 = or 6 — zbmTT, m = 0, 1, ... 


53 


{x, p, t, k) {~x, —p, t + T/2, —k) 


U{x + L/2) = -U{-x + L/2) and E{t + T/2) = -E{t) 
requires 

s = (or 6lp = ±rmv/2, m = 0, 1, ...) and £2=0 


54 


(x,p,t,K) {x + L/2, -p, -t + 2ts, -k) 


E{t + ts) = -E{-t + ts) and U{x) = + i/2) 
where ts = or T/2, requires 
s = and £2 = (or 6* = ±m7r/2, m = 0, 1, ...) 



TABLE I: Symmetry transformations which reverse the current in the system (|2I4|I . 



of temporal period T — 27r/aj. Here the time G i'^iT] 
indicates the switch-on time of the driving force E{t). 
We use L/2tt, {ni/kj^UoY/^ , ^"^^ as the units of dis- 
tance, time and energy, correspondingly [s^]- The sys- 
tem given by Eqs. ([2]|4]) could describe the dynamics of a 
diluted cloud of ultracold atoms placed into optical po- 
tential formed by counter-propaga ting laser beams with 
periodically modulated intensities |19l |. 

The Hamiltonian ([2]) is periodic both in time and 
space. The solution of the eigenvalue problem for 
the corresponding Floquet operator U{t,to), \ip{t + 
to)) = U{t,to)\ip{to)), provides the set of cigenfunctions, 
{lipait))}. The cigenfunctions satisfy the Floquet theo- 
rem, \^c.At)) = e-'"°M*/'*|0a^t)), |</'a,«W) = I0a,«(t + 

T)), and the Bloch theorem The Hilbert space of 

the system is sliced into invariant subspaces, each one of 
which is spanned by the states bearing the same quasi- 
momentum value k S [— tt/L, tt/L]. The quasienergy 
values are confined to the interval [—hTr/T^hn/T], and 
quasienergies of a given index a form a band across 
the K-space, eQ[K]. For a given k and theme t the set 
of cigenstates \ipa.K{t)) forms a complete orthonormal 
basis for the corresponding subspace of the total sys- 
tem's Hilbert space. Any initial state, \ip{x,t = to)), 
can be expanded in the Floquet basis, {"tplt = to)) = 

Jli^L'^'^^a'^a,K{to)\4'a,K{to)), and the time evolution 
of the state is governed by 0, [3 • 

\m)= d'^E^".«(^o)e""°"'*^''iV'a,KW). (5) 

The asymptotic current generated by the system. 



J{to) = limj^oo /to takes on a sim- 

ple form in the Floquet basis [Tsj : 

J{to)= / dn}^ Ca,K{to)Va[K], (6) 
•^-1/2 ^ 

where Va [n] is the the average velocity of the correspond- 
ing Floquet state, Valn] = {{(j)aK{t)\p\(l)aKit)))T, and 
(...)t denotes the averaging over one period of the driv- 
ing. Following the Hcllmann-Feynman theorem, the av- 
erage velocity of an eigenstate is equal to the first deriva- 
tive of the corresponding quasienergy band with respect 
to id€c,[K]/dK 0, H^l- The Hilbert space of 

the system can be imagined as a set of "conveyer belt" 
- eigenstates, each one of which is moving with its own 
velocity, and the overall ratchet current depends on how 
the initial wave packet \il>(t = to)) was distributed among 
the various belts [l^ ■ 



III. SYMMETRIES IN QUASIMOMENTUM 
SPACE 

The transport properties of Floquet states are gov- 
erned by the space-time (a)symmetries of corresponding 
Hamiltonian |18| . There are four types of relevant sym- 
metry transformations for the Hamiltonians of the type 
([2]), see Tab. 1. Two of them. Si and S3, demand the 
space inversion, and the remaining two, i.e., ^2 and ^4, 
involve the explicit time reversal. Each transformation 
changes the sign of the quasimomentum, and maps every 
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Floquet band onto itself, a negative branch onto a posi- 
tive one, and vice versa. Si : £a[K] = The Floquet 
states with k ^ may possess nonzero velocities, Uq[k], 
even in the presence of a symmetry; but in this case we 
have 'L'c[[— k] = — Uq,[k]. Note that only the eigenstates 
from the center of the Brillouin zone, Ua[0], are mapped 
onto itself under the transformations from Tab. 1. The 
symmetry analysis of quantum ratchets performed in pre- 
vious works was restricted on transport properties involv- 
ing the zero-quaismomentum states only [l8| . The main 
result of the analysis was that whenever one of the sym- 
metries listed in Tab. 1 holds, all the Floquet states with 
K = are nontransporting, i.e. Va[0] = 0. Independent 
of a particular shape of the initial cloud and the choice of 
starting time the asymptotic current would be absent 

- for any initial state drawn from the subspace k = 0. 
Although such the initial state can serve as a good ap- 
proximation of a diluted BEC cloud smeared over suffi- 
ciently many wells of the periodic potentials |29i], it can- 
not mimic the experimentally relevant situation with a 
partially delocalized wave packet, ip{x,to)- Such a more 
realistic wave packet occupies a certain region in the n- 
space, and therefore, Floquet states of different quasimo- 
menta k and band indices a, will contribute to the wave- 
packet dynamics according to the corresponding weights, 
CaAto), see Eq. ©. 

We could make the reasonable assumption that at time 
to we start out with the initial cloud in form of a symmet- 
ric wave-packet, 'ip{—K,to) = 'ip{K,to)- As it is the case 
with the dynamics restricted to the subspace of k = 0, the 
evolution of a wave packet depends on the initial starting 
time to, Fig. 1(a). In contrast to its classical counterpart, 
however, where the chaotic dynamics erases the memory 
about the initial conditions of the system preparation ex- 
ponentially fast [l8l.[l9|. a fully coherent quantum ratchet 
maintains this information, which is encoded in the co- 
efficients C(to) forever. All the symmetries from Table 
1, except the first one. Si, involve time transformations 

- time inversion in 5*2, time shift in S3, and a combined 
transformation in S'4 - so that the contributions from the 
eigenstates with opposite quasimomenta, might be differ- 
ent for a given, fixed to, Ca^-nito) 7^ Ca^nito), even when 
one (or all) of the symmetries S2, S3, S4 are obeyed. The 
asymptotic current J(to), sec Eq. may acquire then 
a nonzero value. Only the presence of the symmetry Si 
guarantees the absence of the quantum ratchet current 
for any symmetric chosen initial wave packet and for any 
choice of the starting time to- 

The averaging over to provides the way for the re- 
alization of the symmetries 82,83 and S'4, so that the 
equalities for the averaged contributions, {Ca,^Kito))to = 
{Ca,K.{t(j))toi now hold whenever at least one of these 
symmetries is present. The contributions with opposite 
velocities cancel each other out, and the to-averaged cur- 
rent 

J={J{to))u (7) 
then equals zero, see Fig. 1(b). 




FIG. 1: (Color online) Left panel: (a) Quantum ratchet cur- 
rent as the function of the starting time to- The lines corre- 
spond to the asymptotic current calculated by using Eq. ([6]) 
for e = ■k/2 and Op = 2.2 (solid blue); Op = -7r/2 (dashed 
pink). Triangles (squares) correspond to the current calcu- 
lated by the direct propagation of an initial wave packet up 
to the time lOOT. (b) The averaged current J = {J{to))to at 
the spatial phase difference 6p = n/2 versus the phase differ- 
ence 6 between the two harmonic driving components. Right 
panel: The quasienergy spectrum as a function of the quasi- 
momentum. The color plot shows the K-representation of the 
initial wave packet, Eq. ((8|, with ao — 4n. The other pa- 
rameters are cj = 1, ft = 1, Ei = 0.9, E2 = 0.45, K = 0.9, 
s = 0.7, 9 = n/2. Op = ~n/2. The current is in units of the 
recoil momentum. 



IV. QUALITY MEASURE OF QUANTUM 
RATCHET TRANSPORT 

In order to explore the transport performance of 
Hamiltonian quantum ratchets we use here two indepen- 
dent approaches. Namely, we employ the direct inte- 
gration of the Schrodingcr equation ^ on the interval 
[—100 X L, 100 X L] by using a grid with lO'' points, and, 
independently, use the Floquet approach based on the 
plane- wave expansion [l8| . As our initial wave function 
we used a Gaussian wave-packet, reading, 

^(x, to = 0) = (27ra2)-i/4 ^xp (^-^^ , (8) 

with the variance cro = 47r, Fig. [U right panel. 

The results of both numerical schemes exhibit very 
good agreement, see Fig. [TJa). The average current J{6) 
matches the prediction of the symmetry analysis, reveal- 
ing the property J{-9) = -J{0) = J{0 + n) [l^, see in 
Fig.Hb). 

A desirable performance of a quantum ratchet trans- 
port is when the spreading of the wave packet constitutes 
a slow process when compared to the overall directed 
motion of the packet center of mass. This limit would 
correspond to a highly coherent transport, when the bal- 
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listic motion of the atomic cloud is minimally hampered 
by diffusion. The natural measure of the quality of the 
ratchet transport would then be the quotient of the di- 
rected current and a proper diffusive characteristics. The 
diffusion is naturally quantified by using the dispersion 
of the wave packet, cr^(t) = {x{t)'^) — {x{t))'^ . In the case 
of normal diffusion, dispersion scales like cr^(i) oc t, and 
in the classical limit, it perfectly describes the kinetics of 
over- [H, [2^ and under-damped [13] ratchets. Yet 
the situation with normal diffusion is rarely the case even 
in the classical, dissipation-free limit, where ac-driven 
Hamiltonian systems demonstrate typically superdiffu- 
sive kinetics, (t)_oc ^21, with a parameter-sensitive ex- 
ponent 1 < 7 < 2 [30, [3l| ■ In the coherent quantum limit 
the diffusion is uniformly ballistic in the asymptotic limit 
t > ibai, (^"^{i) , with a parameter-sensitive transient 
time scale ibai [lO]: see the inset in Fig. 2(a). Therefore, 
the Peclet number, Pe = JL/D [23|, which previously 
was successfully applied for the calibration of classical 
ratchets, presents now an inappropriate measure in the 
quantum context. We next introduce here its quantum 
analog, named quality of ratchet transport^ 

Q = Q{to) = J''{to)/DMto) , (9) 

where Dcs = limt_j,oo (T^(i)/i^. This so introduced "qual- 
ity" quantifier, as it is the case for the classical predeces- 
sor, Pe, is dimensionless. 

Among the several parameters of the system dynamics 
it is the driving frequency cj which presents a most suit- 
able control parameter. This driving frequency w allows 
one to tune the ratchet into resonant regimes [l^- The 
frequency of the driving can conveniently be varied in ex- 
periments [1^ . The dependence of the effective diffusion 
coefficient Doff on the driving frequency, see Fig. 2(a), 
does not keep abreast with the corresponding current de- 
pendence, see Fig. 2(b) (thin line), therefore producing 
regions of different transport quality. Note that although 
the asymptotic absolute current values | J| are identical 
at the driving values uj ~ 1 and uj = 2.68, the quality 
of the transport in the first point is almost three times 
better than in the second. 

In order to gain insight into mechanisms of highly- 
coherent quantum transport, we resort to the analysis 
of the system Floquet spectrum. First we sorted out all 
the Floquet bands according to their kinetic energy, and 
then ordered them into an ascending order. This is a 
reasonable procedure, since the states with kinetic ener- 
gies that are much larger than the characteristic energy 
KUoi would not contribute tangibly to the initially pre- 
pared wave packet. Finally we left with a small number 
of relevant Floquet bands of low energy, see Figs. 3(a,b). 

The relative slopes of the energy bands determine not 
only the velocity of Floquet eigenstates, but also the dif- 
fusionproperties of the system. By using the results of 
Rcf. we can rewrite the expression for the effective 
diffusion coefficient in the following form (the detailed 




UJ 



FIG. 2; (Color online) (a) Dependence of the effective diffu- 
sion coefficient Dca ~ Dca(to — 0) on the driving frequency uj 
for the initial wave packets, Eq. [51 of the dispersion ao = in 
(thin black line) and ao = 407r (thick red line). The inset 
depicts the time evolution of the dispersion (t) of the wave 
packet with ctq = 47r for uj — 1. The dashed line corresponds 
to the asymptotic limit cr^ (t) oc t'^ . (b) The quantum ratchet 
current J{to — 0) and (c) the quality of the coherently di- 
rected quantum transport Q{to = 0) as functions of the driv- 
ing frequency cj. The two circled values at a; = 1 and at 
(jj — 2.68 are discussed further as functions of quasienergy- 
quasimomentum characteristics in Figure 3. The diffusion 
coefficient and the current are in units of the square of the 
recoil momentum and the recoil momentum, correspondingly. 



derivation is given in the Appendix), 

1 /"/"^/^ -s-^ 
^cff(^o)=2 // dKidK222\'^^,-i(io)\^\'^'P,^2ito)\'' 

JJ-l/2 

^ / dea[Ki] ^ dep[K2\\'^ 
\ dni dK2 J 

(10) 

With the Eqs. ^ and (fTO| in hands, we can discuss 
two extreme limits. Consider first the ultimate case with 
only two bands, i and j, that are populated initially. If 
the bands are parallel at the vicinity of k = 0, both 
have non-zero slopes and share the initial cloud equally, 
Cj(^o)P = |C'i(^o)l^ ~ 1/2; then the corresponding effec- 
tive diffusion coefficient would be equal to zero. Accord- 
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-0.2 -0.1 0.1 0.2 -0.2 -0.1 0.1 0.2 



K, K 

FIG. 3: (Color online) The quasienergy bands Sa as a func- 
tion of quasimomentum k. (a) The case of high-quality quan- 
tum ratchet transport realized at the circled value co — 1 
in Fig. 2(b). (b) The case of a low-quality quantum ratchet 
transport as realized at the circled value a; = 2.68 on Fig. 2(b). 
The depicted "widths" relate the color coding [from dark 
(maximal) to bright (minimal)] of a band encode to the corre- 
sponding weig ht, ICc^to = 0)p. The weights were obtained 
by projecting the initial Gaussian wave packet ((8]) onto the 
corresponding Floquet state 4>a,K{to ~ 0). 



ingly, the quality Q would be infinite in this case. This 
regime is approached at the frequency uj = I (point of a 
high-quality transport in Fig. 2(c)), see Fig. 3(a). 

The situation at the frequency lu = 2.68 is contrary to 
the discussed one, see Fig. 3(b). There we also find only 
two mostly populated bands which exhibit an avoided 
crossing (AC). The bands have different slopes, and ac- 
cording to the last factor in Eq. ([TU)) . this implies strong 
diffusion. Both bands are almost horizontal within the 
AC region, so that they yield almost no current. Most 
of the contribution to the overall current stems from the 
parts of the bands, located away from the AC region. 
Although the bands have a large slope, thus producing a 
directed current which matches that of the regime at the 
point UJ = 1, the presence of strong diffusion decreases 
the transport quality by factor of three (point of a low- 
quality transport in Fig. 2(c)). 



V. CONCLUSION 

In summary, we have investigated the coherent quan- 
tum transport of a quantum wave packet in a flashing 
periodic potential. The regime of high-quality ratchet 
transport corresponds to the ballistic center-of-mass mo- 
tion minimally obstructed by quantum diffusion. The 
quality is a tunable characteristics, and it is indeed pos- 
sible to optimize the quantum ratchet performance by 
tailoring the parameters of the driving flashing poten- 
tial. The transport quality Q is also an experimentally 
relevant characteristics: the ratchet current can be mea- 



sured by using the time-of-flight technique [l^ , while the 
dispersion of BEC cloud can be estimated directly from 
absorption images [3^ . 



The quality of the transport depends on the initial 
time to, since the band populations, |CQ^„(io)P, are func- 
tions of starting time to. Therefore, different choices 
of the starting time to may lead to regimes of differ- 
ent transport quality even when all the driving param- 
eters of the system Hamiltonian are fixed. The depen- 
dence on the dispersion of the initial wave packet, (Tq 
([5]), is even more stringent. Namely, the localization 
in the x-space is inversely proportional to the localiza- 
tion in K-space. Therefore, a very delocalized initial 
wave packet, smeared over many periods of the opti- 
cal potential, would be strongly localized at the point 
K = in the quasimomentum space. The asymptotic 
transport dynamics of the wave packet can be evaluated 
from the behavior of the quasienergy bands at the cen- 
ter of the Brillouin zone, by using Eqs. (|6]) and (|10p . 
The localization in K-space, thereby, makes the effect 
sharper and typically leads to the enhancement of the 
quality of ratchet transport, see Fig. 2(c). On the con- 
trary, the ratchet effect would be faint for an initial wave 
packet in the form of a very narrow peak in x-space. 
Strong localization of the cloud would correspond to an 
almost uniform distribution over the Brillouin zone, and 
since the average band velocities are all equal to zero, 
J^l^j^ Va[K]dk = ea[—TT/L] — CaliT / L] = 0, such the initial 
distribution will yield almost imperceptible current. At 
the same time the spreading rate would be huge. Para- 
doxically enough, a very narrow wave packet will lead to 
quantum ratchet transport of a very poor quality [33| . 



Our results can be considered as a first step towards an 
optimization scheme of cold atom engines. The ultimate 
aim of the research in this direction is an operational 
recipe, a protocol, which would allow one to maximize 
the quality of a quantum ratchet transport without 
going deeply into the analysis of the particular system 
of interest. One of the intriguing subdirections is the 
role of disorder in the functioning of quantum ratchets. 
Namely, a following question is of special interest: Can 
a pinch of disorder in the underlying optical potential 
improve the performance of quantum ratchets? In 
other words, can the disorder- induced localization [s^l 
slow down the diffusive spreading without affecting the 
center-of-mass velocity? 
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Appendix A: EFFECTIVE DIFFUSION 
COEFFICIENT 



We start out from the expression for the mean square 
dispersion given in the Appendix B of Ref. , 



it f ' dKV|a,(t 

V ^-1/2 a 



o;i'%i^l +0(t). (Al) 



Here and in the following we drop the dependence on the 
initial starting time to for the results on l.h.s. of equa- 
tions. The corresponding effective diffusion coefficient is 
given by, 

We consider the integral with respect to k as the limit of 
the summation, 



/1/2 n 
dn — >■ lim } An 
1 /, Ak->0 ^-^ 



(A3) 



with the condition tiAk = 1. Next, the expression (jA2[) 
is transformed into a summation with respect to indices 
K and a. Wc reorder all the terms from 1 to A^, and 
recast Eq. (jA2p in the following form. 



N / N 

„2 



Dcff = ^ aix\ - ^ aiXi 

i=\ \i=l / 



(A4) 



where = AK|CQ,K(to)P and Xi = d£a\K\ldK. 

Due to the norm preservation, we obtain X^iLi — 
1. We further simplify Eq. (|A4[) . and by using a simple 
algebra, arrive at 



N 



N 



N 



DcS = ^ aiXi - ^ aiajXi 

N 

= '^^ai{l ~ ai)x1 ~ 2'^^aiajXiXj (A5) 

1=1 i=l 
j>i 

1 ^ 

= - a^aJ{x^ -Xj)"^. (A6) 

We next return to the summation with respect to the 
indices k and a and take the limit An — > 0. The integral 
form is recovered and therefore we end up with Eq. (|10p . 
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